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Abstract 

It is well known that the Hubbard model on a line graph has a flat band and ferromag- 
netic ground states in a certain density range. We show that for a Hubbard model on a 
line graph of a planar bipartite graph the ferromagnetic ground state is stable if one adds 
a special contribution to the kinetic energy which lifts the degeneracy of the lowest single 
particle state. Stability holds for sufficiently strong repulsion U . The model has extended 
single particle eigenstates, no degeneracy, and no band gap. It is therefore a good candidate 
for metallic ferromagnetism. 

PACS-numbers: 75.10.Lp, 71.10.Fd 



1 Introduction and Results 

Ferromagnetism in itinerant electron systems is an old and not yet fully understood problem in 
statistical physics. Kanamori [l] was the first to point out that the Hubbard model is a good 
starting point for the investigation of ferromagnetism. It contains the two essential ingredients 
one needs to understand itinerant ferromagnetism: A single particle hopping term describing the 
motion of the electrons on a lattice and a purely local, repulsive interaction U between them. 

The first major result towards ferromagnetism in the Hubbard model was the so called Theorem 
of Nagaoka It states that for an infinite local repulsion and very close to a half filled band 

(one electron in addition to half filling) the ground state of the Hubbard model is fully polarized. 
The most general proof of this theorem was presented by Tasaki [Hj. It uses a special basis for 
the multi-particle states in the Hilbert space and the Perron-Frobenius Theorem for matrices. 

But the Theorem of Nagaoka describes the situation close to the singular point where the repulsive 
interaction is infinite and the band is half filled. In that case, each lattice site is occupied by one 
electron and one has a complete spin degeneracy, a pure paramagnet. Therefore, the Nagaoka 
Theorem is not a good starting point for a theory of ferromagnetism in the Hubbard model. 

In 1989, Lieb jSj proved a result on the Hubbard model on bipartite lattices which yields a 
ferrimagnetic, not fully polarized ground state if the number of lattice sites in the two sub- 
lattices of the bipartite lattice is different and the difference is an extensive quantity. One can 
show that in this case the model has a purely dispersion-less or fiat band and that there are other 
bands symmetrically above and below the fiat band. Half filling means that the Fermi energy 
lies in the fiat band, which means that the density of states at the Fermi level pp is infinite. 
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This result is in some sense complementary to the Nagaoka Theorem. The old criterion of Stoner 
for the occurrence of ferromagnetism, which is a simple mean-field result, states that ppU > 1 for 
ferromagnetism to appear. Here U is the repulsive on-site interaction between electrons and pp 
is the density of states per site at the Fermi level. This criterion over-estimates the occurrence of 
ferromagnetism. One knows examples of special lattices where even for pp = oo and U = oo the 
ground state is not ferromagnetic. Nevertheless, an infinite value of pp may be a good starting 
point for the study of ferromagnetism in the Hubbard model. 

In 1991, the first example for a Hubbard model with pp = oo was published which has a unique, 
fully polarized ferromagnetic ground state, the Hubbard model on the Kagome lattice [3 EI- This 
kind of ferromagnetism is called flat band ferromagnetism. It was shown to occur for a large class 
of lattices which are line graphs. Tasaki [niE3 proved the existence of flat band ferromagnetism 
on another class of decorated lattices. Later on, a general criterion was derived which yields 
a necessary and sufficient condition for the uniqueness of the ferromagnetic ground state for a 
general situation where the hopping term in the Hubbard model has degenerate single-particle 
ground states This result applies to any lattice with a lowest flat band, it applies even 

in a single band case with a partially flat band, but one needs a special lattice topology which 
is expressed as a criterion on the single-particle density matrix. Sekizawa recently proved 
an extension of Tasaki's result for a general class of decorated lattices, which fulfill the criterion 
derived in [T2|. 

Line graphs of planar bipartite graphs have special properties. The reason will become clear 
below. Using these special properties, we have been able to give a complete description of the 
ground states of the Hubbard model on this class of graphs in a certain density range pi]. Here 
we make use of these special properties as well. 

Often, it was argued that the flatness of the band is artificial and that flat bands do not occur 
in nature. Another point is that in a flat band the Wannier states are single particle eigenstates 
of the Hamiltonian. The Fermi energy lies in a region where model has localized single-particle 
eigenstates. Although, due to the degeneracy, the distinction between localized and extended 
states is irrelevant, is is clear that arbitrary small local perturbations (point defects, line defects, 
disorder) favour localization. This means that in a realistic situation one probably has a ferro- 
magnetic insulator and not a ferromagnetic metal. This argument does not apply to the case of 
a partially flat band, but even than the extensive degeneracy of the band is artifical. 

Therefore, already in the beginning of the discussion of flat band ferromagnetism, Tasaki tried 
to show that the flat band ferromagnetism is stable if one introduces a small dispersion to the 
flat band and if the repulsion is sufficiently large. A first rigorous result in this direction was 
published in 1995 for a special class of decorated lattices where there is a gap between the 
flat band and the other bands ^3 El- One year before, Kusakabe and Aoki published some 
numerical results that support the stability of flat band ferromagnetism on small lattices [T7jll8|. 
Finally, in 2003, Tanaka and Ueda showed that even on the Kagome lattice, which has no gap 
between the flat band and the next band, flat band ferromagnetism is stable if one adds a special 
hopping term to the Hamiltonian that yields a small dispersion in the formerly flat band |19j . 

All these examples indicate that at least for a large class of lattices flat band ferromagnetism 
may be stable: If one adds additional hopping terms to the Hamiltonian so that the flat band 
gets a dispersion, the ferromagnetic ground state remains the unique ground state for sufficiently 
large U. The present paper is a step in this direction. 

Recently, Tamura et al proposed to study flat-band ferromagnetism experimentally using quan- 
tum dot super-lattices (201 l^- They performed calculations for a Kagome lattice within a 
tight-binding approximation for the quantum dot super-lattices, which includes a next nearest 
neighbour hopping term. The model was solved using exact diagonalization of small lattices. 
Although the next nearest neighbour hopping destroys the flat band, the result shows that flat 
band ferromagnetism can be observable in dot arrays. 

Let us now summarize the results of the present paper: 
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Our results applies to a large class of graphs, line-graphs of planar, bipartite, two-connected 
graphs. It is the same class of graphs as in This class contains many two-dimensional 

lattices, but we do not need translational invariance. One property of the nearest-neighbour 
hopping on such graphs is that the lowest single-particle eigenstate is highly degenerate. The 
ground state of the full Hamiltonian for the case where the number of electrons equals the 
number of degenerate single-particle ground states is ferromagnetic and unique up to the usual 
spin degeneracy. We add a contribution to the kinetic energy which lifts the degeneracy. Then 
we are able to show that the ferromagnetic ground state remains the unique ground state of the 
system if the repulsive interaction U and the nearest-neighbour hopping t are sufficiently large. 

The result of the present paper is a generalization of the result in ^H] for the Kagome lattice to a 
large class of lattices. In fact, the experimental realizations discussed in jSHlI^ are line graphs of 
finite planar bipartite graphs. Furthermore, the additional hopping term contains many (local) 
parameters. And the result may be a good starting point for more general results in the future. 
For the Kagome lattice, our proof is a simplification of the proof in jl9j . 

The main difference between our result and the results obatined by Tasaki is that in our 

case there is no band gap between the lowest band and the higher bands. In addition, at least in 
the translational invariant case, the single particle eigenstates on our lattices are not localized. 
Localized single particle eigenstates or a gap in the single particle spectrum are criteria for a 
Mott or Anderson type of insulator. Therefore our class of graphs contains lattices which are 
good candidates for metallic ferromagnetism. 

The rest of the paper is organized as follows. In the next section we introduce the notation, some 
definitions, the conditions we need, and we state the main result. We also give some comments 
on that result. The third section of the paper contains the proof. 

2 Definitions and Results 

As stated in the introduction, we investigate the Hubbard model on line graphs of planar bipartite 
graphs. Let us therefore introduce some definitions and properties concerning line graphs of 
planar bipartite graphs. We use some notions of standard graph theory which may be found in 
any text book on graph theory, see e.g. [221 ■ ^^e the notation in 

Let G be a bipartite planar 2-connected graph. V{G) is the set of vertices, E{G) the set of edges 
of G. Each edge e S E{G) connects two vertices x,y £ V{G) and may be denoted as e = {x,y}. 
As in 23] we assume that G is 2-connected which means that no edge exists so that the graph 
decays into two unconnected subgraphs if this edge is deleted from G. A graph is bipartite if 
the vertex set ^ is a union of two disjoint vertex sets Vi and V2 so that each edge joins a vertex 
from Vi to a vertex of V2. In a bipartite graph, the length of each cycle is even. 

A walk of length n — 1 is a sequence c = (xi, ei,X2, 62, ... , e„_i, A path is a self avoiding 
walk, a cycle is a closed, self avoiding walk. 

A planar graph can be represented in the plane as a set of point and lines between them so that 
there are no lines intersecting each other. The representation of the planar graph in the plane is 
called a plane graph. The representation is not necessarily unique. In the following we assume 
that the planar graph G is a plane graph. The results do not depend on the representation in 
the plane. 

A plane graph divides the plane into a set of connected components, called faces. Most of the 
faces are bounded. Each plane graph has exactly one unbounded face. Let F{G) be the set 
of bounded faces of G. Due to the famous Theorem of Euler, the number of bounded faces is 
|F(G)| = \E{G)\ - \V{G)\ + 1. We assume \F{G)\ > 1. If there is no or only one bounded face, 
the results we want to show are trivial and uninteresting. The boundary of a face / is a cycle. 

Let L{G) be the line graph of G. The line graph of G is defined as follows: The vertex set of 
L(G) is the set of edges E{G) of the original graph. The vertices e,e' E E(G) are connected if 



3 



|e n e'l = 1, i.e. if the two edges of the original graph have a vertex in common. Note that a 
hne graph of a plane graph is not necessarily a plane graph. The Kagome lattice we mentioned 
above is the line graph of the hexagonal or honey-comb lattice. It is a plane graph. But the line 
graph of the square lattice is not planar. 

Let B = {bxe)xeV{G), eeE{G) tie the vertex-edge-incidence-matrix, bxe = 1 if the vertex x belongs 
to the edge e, bxe = otherwise. 

Let B = (bfe) f^F{G),e£E{G) t>e the face-edge-incidence-matrix, bj^, = 1 if the edge e belongs to 
the boundary of /, 6je = otherwise. 

The number of edges of the boundary of / is = J2eeE{G) ^/e- Since the graph is bipartite, nj 
is even. 

Let C = {cxf)xeV{G), f€F{G) be the vertex-face-incidence-matrix. Cxf = 1 if the vertex x belongs 
to the boundary of /, Cxf = otherwise. One has 

Cx/ = ^ ^/e^^e (1) 

ee-B(G) 

We let dx = J2f£F{G) '^xf and D = diag{dx)- dx is the number of faces touching the vertex x. 

Each face / G F{G) may be oriented. We choose a clockwise orientation. Each e may be oriented. 
Since G is bipartite, we choose the orientation of e such that it points from Vi{G) to V2{G). We 
now introduce S = {s fe) f^F(G),e&E{G) where s/e = 1 if e is part of the boundary of / and points 
into the direction of /, sje = — 1 if e is part of the boundary of / and points into the opposite 
direction of /, s/e = otherwise. 

A standard result in graph theory is BS* = 0. Furthermore, the dimension of the kernel of B is 
|^(G)| - \V{G)\ + 1 = \F{G)\. The columns of 5* form a basis of the kernel of B. 

We now consider a Hubbard model on L{G) with the Hamiltonian 

H = ^ tee'cLceV + Y UeUe^Uei (2) 
e,e'G£(G),cr ee-B(G) 

where 

tee' = t'^bxebxe' — X! ~^}eSfe'- (3) 

We let f > 0, s/ > and Ue > 0. c\„ and Cea are the usual creation and annihilation operators 
on the vertices of L{G). Uea = c^^Cea is the number operator for electrons on e with the spin a. 
N = X^eo- ^eo- is the number of electrons. The first part of the Hamiltonian describes the hopping 
of electrons on L{G), the second part is the local repulsive electron-electron interaction. 

The first part of tee' is the usual nearest neighbour hopping on L(G) plus an additional term 2t5ee' ■ 
The lowest eigenvalue of ^a:e^xe' is zero. It's multiplicity is 1-^(0)1 = — -|- 1. The 

second term lifts the degeneracy of this eigenvalue and does not change the other eigenvalues. 
Since the second part of the matrix {tee') is negative semi-definite, all eigenvalues which are 
for = become non-positive. For s/ = this is the model discussed in |14j. 

The second part of the kinetic energy of the Hamiltonian can be written as i?i = — X]/o- ^/^/o-'^/o- 
where = 'n^^^'^Yl,e^fec\a- Since the corresponding single particle states are not orthogonal, 
^//' = = {iT'fnf y^/'^J^eSfeSf'e is not diagonal. 

Let Sp = diag(sj) and S = (s fg) f,g£F{G)- Then single particle spectrum of Hi is given by the 
spectrum of the matrix —-/SfS^'^Sf. 
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The Hubbard model has the usual SU (2) spin symmetry. It commutes with the spin operators 

= ^ (™eT - nei) (4) 

eeE(G) 

e(^E(G) 

S- = 4iCeT (6) 

e(^E(G) 

The eigenstates oi H are eigenstates of S"^ = ^(S'+S'_ + S'_5+) + S'2 with eigenvalues S'(S' + 1) and 
Sz- Due to the SU{2) symmetry, an eigenstate of H with a spin S is {2S + l)-fold degenerate, 
5*^ takes values between —S and S for the degenerate states. 

Let us now state our main result: 

Theorem: For sufficiently large t > tc and U > Uc (both independent of the number of lattice 
sites and finite in the thermodynamic limit) and for A'' = |-F(G)| electrons, with G a 2-connected 
planar graph, the ground state of the Hubbard Hamiltonian Q on the line graph L{G) is fully 
polarized, S = ^N, and unique up to the usual (25 + l)-fold spin degeneracy. In the subspace 
Sz = S the unique ground state is H/ "^/flO) where |0) denotes the vacuum. 
Before we start to prove this theorem, let us make some remarks: 

(i) The ground state in the theorem is a fully polarized state. It's energy does not depend on U. 
This state is an eigenstate of the Hamiltonian for any value of U. But it is clear that it is not 
the ground state for U = 0. For ?7 = we obtain the ground state by filling the single particle 
eigenstates of Hi successively with electrons with spins up and down. An exception occurs if 
a single particle eigenstate is degenerate, in that case Hund's rule applies. But generically, the 
degeneracy will not be of the order \F{G)\. 

(ii) The proof of the theorem makes use of the fact that the columns of S* form a basis of the 
kernel of B. In ^1] we made use of this fact to completely describe all ground states of the 
Hamiltonian with Sf = for N < \F{G)\. Now we use this property to prove a result for the 
ground state of H for Sf > and = \F{G)\. One might expect that it should be possible to 
extend this result to < \F{G)\. But this is much harder. Even if s/ = s for all / G F{G), a 
simple generalization of the theorem is not possible. The reason is that we use H > —J2f for 

sufficiently large U and t together with the fact that the ferromagnetic state ri/^/flO) has the 
energy —J2f^f- For < |-F(G)| a generalization is difficult, since the proof uses the locality of 
the states dj^. But the single particle ground states of Hi are not local and therefore it is not 
possible to express the ferromagnetic eigenstates of H with lowest energy by local operators. 

(iii) Although the proof makes use of the locality of the operators dj.^, the single particle eigen- 
states of H are generically not localized. Esp. this is not the case if G is a translational invariant 
lattice and sj = s. In that case the system is a good candidate for a ferromagnetic metal: There 
are no artifical degeneracies, the single particle eigenstates are extended states, not localized 
states, and the Fermi energy does not lie in a band gap. 

(iv) The proof makes use of the fact that G is a plane graph. For arbitrary line graphs, the 
degeneracy of the single particle ground state is — |y| + 1 if G is bipartite, as for the case of a 
plane graph. One can even construct a basis using small cycles on G, similar to the boundaries 
of the faces of the plane graph. But the number small cycles is much larger than the number 
of single particle ground states. As a consequence, the bounds we shall use in the proof do not 
hold. Therefore, it is not possible to easily generalize the result to three dimensional lattices. 

3 Proof 

The Hamiltonian H can be written as a sum of local Hamiltonians H = J2 /gF(g) where 
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(7) 



Hf,l = -zf- Sfecl^Sfe'Ce'a (8) 

e,e'eE(G),CT 

^/,2 = ^ 13 ^ X] ^a;e4a^a:e'Ce'(T (9) 
xeV{G) ^ e,e'&E{G),a 

^/.3 = H fexeC^e^e+ne- (10) 

Let us first discuss Hf. One has > —2sf. We will show that for sufficiently large t and Ue 

one obtains iJ/ > —Sf. 

In a first step, we discuss Hj for Ue = oo and t = oo. i7/-_2 > and t = oo means that the 
single Hilbert space is restricted to single particle states with Hf^icp = 0. Tlie single particle 
Hamiltonian Hf 2 is equivalent to the matrix 

Tf,,=tB'CfD-^CfB (11) 

where Cf = diag(ca;/) and D = diag((ia;) has been defined above. The kernel ofTf^2 is the kernel 
of CfB. 

In the following discussion we restrict ourselves to the subset of edges of E{G) for which at 
least one matrix element of Tf^2 does not vanish. This subset consists of the edges of / and 
the edges connected to /. One can show that the dimension of the rank of CjB is j/j — 1 if 
no other edges are connected to /, |/| otherwise. The first case has been excluded above since 
G must be 2-connected. This means that the dimension of the kernel of CfB is given by the 
number of edges connected to / but not an element of /. We call them outer edges of /. Their 
number is given by rj = J2xef ^-ef ~ ^l/l- Therefore we need to find rj linear independent states 
that form a basis of the kernel of CfB. The first state is ^o(e) = Sfg. The other basis states 
are constructed as follows: We choose eg to be one the outer edges of /. The other edges are 
numbered as ei, . . . , Cr^-i. We now let pi be the path that starts at cq, ends at ej and runs 
clockwise around /. The path is oriented, it starts at cq- If is connected to the same vertex 
in / as eo, it has length 2, i.e. it contains no edge of /. We let ^i(e) = 1 if e belongs to pi and 
has the same orientation as pf, (f)i{e) = — 1 if e belongs to Pi and has the opposite orientation; 
(l)i{e) = if e does not belong to pi. Obviously, we have constructed ry states. It is easy to see 
that eacli of these states belongs to the kernel of CfB. Furthermore, there is at least one edge 
in / for which all (pi vanish except (po. Since ^j, i > is the only state that does not vanish on 
Cj, the states (pi are linearly independent. Thus, they form a basis of the kernel of Tf^2- 

Now, we introduce 

4.=j2^i(^yi (12) 

e 

One has Oq^ = d^f^. Any multi-particle state which is an eigenstate of Hf^2 with eigenvalue 
can be written as 

where V'o.i,! are build by linear combinations of products of a suitable number of aj^, z > 
acting on the vacuum. 

Now we use that U = co. This condition means that ip has no double occupancy. We noted 
above that there exists one edge e for which all (pi vanish except Obviously, the last term in 
ip has a double occupancy on this edge. Therefore, for C/ = co, ipi must vanish. 

Next, we want that ip is a lowest energy eigenstate of Since if/^i^o,±,i = 0, we obtain 

Hf,iip = —Sfip if and only if ipo vanishes. Thus, for ?7 = oo and t = <x> one has Hf > —Sf and 
the ground state of Hf is of the form 

^ = 4.V+ + ^iV-- (14) 
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By continuity, this result holds for sufficiently large but finite U and t as well, i.e. for U > Ucf 
and t > tcf- Ucf and t^f are finite numbers and do not tend to infinity when G becomes large, as 
long as / remains finite. If e.g. G is some part of an infinite lattice, Ucf and tcf do not change 
when one takes the thermodynamic limit. 

As a consequence of the above result one has H > — J2f£F(G) 

Let us now consider the fully polarized state 

/eF(G) 

The operator dg^j = '}2f{S^^)gfdfcr is the dual operator to dj.^; we have [Jpo-, dj-^] = 5fg. We 
now obtain 





(16) 


ZT'^fgd^f^dg^M 


(17) 


yuf 


(18) 




(19) 



/ 



This shows that for sufficiently large values of t and Ue the fully polarized state il)o+ is a ground 
state of H . The uniqueness of •i/'o+ can be shown quite easily. In fact, we proved the uniqueness 
in 25 for = 0. This proof can be carried over without any changes and it is not necessary to 
repeat it here. 
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